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Abstract 

In this paper we present the design of a straight drift- 
quadrupole system which can transport certain beam el- 
lipses (apochromatic beam ellipses) without influence of 
the second and of the third order chromatic and geometric 
aberrations of the beamline transfer map. 

INTRODUCTION 

A straight drift-quadrupole system can not be designed 
in such a way that a particle transport through it will not 
depend on the difference in particle energies and this de- 
pendence can not be removed even in first order with re- 
spect to the energy deviations. Nevertheless, the situation 
will change if instead of comparing the dynamics of in- 
dividual particles one will compare the results of track- 
ing of monoenergetic particle ensembles through the sys- 
tem or will look at chromatic distortions of the betatron 
functions appearing after their transport through the sys- 
tem. From this point of view, as it was proven in [ 1 1, for 
every drift-quadrupole system there exists an unique set of 
Twiss parameters (apochromatic Twiss parameters), which 
will be transported through that system without first order 
chromatic distortions. In this paper we continue the study 
of the apochromatic properties of straight drift-quadrupole 
systems and present the design of a beamline which can 
transport apochromatic beam ellipses without influence of 
the second and of the third order chromatic and geometric 
aberrations of the beamline transfer map. 

VARIABLES, MAPS AND APOCHROMATS 

We consider the beam dynamics in a magnetostatic sys- 
tem and, as usual, take the path length along the reference 
orbit r to be the independent variable. We use a com- 
plete set of symplectic variables z — (x,p x ,y,p y ,a,e) T 
as particle coordinates. Here x, y measure the transverse 
displacements from the ideal orbit and p x , p y are trans- 
verse canonical monenta scaled with the constant kinetic 
momentum of the reference particle po. The variables a 
and e which describe longitudinal dynamics are 

a = cPo {t -t), e = {£ - fi ) / (/3 2 fib), (1) 

where £o, /?o and to = to(r) are the energy of the refer- 
ence particle, its velocity in terms of the speed of light c 
and its arrival time at a certain position t, respectively. 

We represent particle transport from the location r = 
to the location r = I by a symplectic map Ai and assume 
that the point z — is the fixed point and that the map 
A4 can be Taylor expanded in its neighborhood. Addition- 
ally we assume that the transverse motion is dispersion free 
(always true for the drift-quadrupole systems) and uncou- 
pled in linear approximation, which is a restriction on the 
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form of the 6x6 symplectic matrix M which represents 
the linear part of the map M. 

Let go be some function of the variables z given at the 
system entrance. Then its image gi at the system exit under 
the action of the map M is given by the following relation 

Vz 9l (z) = go(M-\z)), (2) 

which symbolically we will write as gi = : M. -r 1 go- 
Let us consider some Courant-Snyder quadratic forms 

f Ix{t) = 1x{t)x 2 + 2a x {r)xp x + /3 x (r)pl 

< (3) 

{ L y( T ) = ly( T )v 2 +2a y {T)yp y + P y {T)p 2 y 

We say that the map Ai is an m-order (m > 2) apochro- 
mat with respect to the incoming Courant-Snyder quadratic 
forms I x (0) and I y (0) if 

: M r 1 4,3,(0) - : M T 1 I x , v {0) = O {\z\ m+2 ) , (4) 

i.e. if the chromatic and geometric distortions to the shapes 
of the ellipses I x ,y(l) after the system passage come only 
from nonlinear map aberrations which are of the order 
m + 1 and higher. We call the Twiss parameters that enter 
the Courant-Snyder quadratic forms satisfying (3) apochro- 
matic Twiss parameters. 

Up to any predefined order m the aberrations of the map 
M. can be represented through a Lie factorization as 

: M : = n exp(: F m+1 • ... • F 3 :) : M :. (5) 

where each of the functions is a homogeneous poly- 
nomial of degree k in the variables z and the symbol = m 
denotes equality up to order m (inclusive) when maps on 
both sides of Q are applied to the phase space vector z. 

Using the representation ((5} we can state that the mag- 
netostatic system is an m-order apochromat with respect 
to I x (0) and I y (0) if, and only if, all homogeneous poly- 
nomials in the Lie exponential representation (0 of 
the system transfer map can be expressed as functions of 
I x {0), I y (0) and e only. 

NONLINEAR ABERRATIONS OF 
STRAIGHT DRIFT-QUADRUPOLE CELL 

The straight drift-quadrupole cell is a magnetostatic sys- 
tem which is symmetric about both, the horizontal mid- 
plane y — and the vertical midplane x — 0. It means that 
in the Lie exponential representation of the cell map 

: M c : = rn exp(: F c m+l ■ ... ■ .11 :} : (6) 

the polynomials J 7 ^ do not depend on the variable a and 
are even functions of the variables {y,p y ) and of the vari- 
ables (x,p x ). If the drift-quadrupole cell is not a pure 
drift space, then the structure of the polynomial F^, which 
is responsible for the second-order map aberrations, can 
be further clarified using the concept of the apochromatic 
Twiss parameters. Let /3" (r), a%^{r) and 7° >2/ (r) be 



the cell apochromatic Twiss parameters and let J° (r) be 
the corresponding Courant- Snyder quadratic forms. Then 
the polynomial T 3 can be written as follows 

H = -{e/ 2) ■ ■ I a x (0) + UPy) ■ ly(0) 

-^•(e/7o) 2 ), (7) 

where l c is the cell length, 70 is the Lorentz factor of the 
reference particle, and £ x (fl%) and £, y (f3y ) are the cell 
chromaticities calculated for the apochromatic Twiss pa- 
rameters (see more details in |Q][3]). 

CONCEPTUAL SOLUTION FOR 
THIRD-ORDER APOCHROMAT 

Let us consider a system constructed by a repetition of 
n identical drift-quadmpole cells and let us assume that the 
horizontal and vertical transfer matrices of the n-cell sys- 
tem are equal to the plus or minus identity matrices, but 
the horizontal and vertical cell focusing matrices are not 
equal to the plus or minus identity matrices. Then, as it was 
shown in 0, the shape distorting effects of the second- 
order aberrations of the cell map on the transport of the cell 
periodic Courant-Snyder quadratic forms are washed out at 
the exit of the 71-cell system by averaging. In other words, 
it means that under the assumptions made the ?i-cell sys- 
tem is a second-order apochromat with respect to the cell 
periodic Courant-Snyder quadratic forms. And, as we will 
show below, the idea of the usage of averaging provided 
by repetition of n identical cells for automatic cancellation 
of undesired effects can also be applied to the design of a 
third-order apochromat if the cell periodic Twiss parame- 
ters coincide with the cell apochromatic Twiss parameters. 

So, let us start with the assumption that the cell transfer 
matrix allows periodic beam transport and that the cell pe- 
riodic Twiss parameters are, at the same time, the cell apoc- 
hromatic Twiss parameters. Then the map of the repetitive 
71-cell system can be expressed as follows 

: M nc := 4 exp(: S 3 + S 4 + S 5 + S 5 :) : M c " :, (8) 

where the homogeneous polynomials (aberration func- 
tions) S3 , S4, 1S5 , 1S5 are given by the formulas 

S 3 (z)=n-^(z), (9) 

S 4 {z)=n-n{Ti{z)), S 5 (z)=n-TZ(^(z)), (10) 

S 5 (z) = (l/2){j^(z),W n (z)}, (11) 
the binary operation {*,*} is the Poisson bracket, 

1 n— 1 

*(/(*)) = - E /( M c m *)> < 12 ) 

771=0 

[n/2] 

Wn(z) = (n + l-2m) 

■ (Tl{M r c l - m z) - Fl(M™- x z)), (13) 

and the symbol [k] denotes the biggest integer which is 
smaller or equal to k. 



Our second assumption is that the periodic cell phase 
advances /i£ and fi y satisfy 

H° xy = 2nq Xt y I n (mod27r) (14) 

for some integer or half integer q x and q y which are smaller 
than n and are such that the resonances 

2/4 iy , 4/4^, 2/4 ±2/4 (15) 

are avoided. From this assumption it follows that the func- 
tions 54 and 55 (as a result of averaging provided by the 
operator TZ) can be expressed as functions of the variables 
J" (0) and e only, and thus the n-cell beamline becomes 
a third-order apochromat with respect to the cell periodic 
Courant-Snyder quadratic forms I£ (0). 

Note that the minimal number of cells required for the 
construction of the third-order apochromatic beamline fol- 
lowing the above recipe is five, and that the first distortions 
to the shapes of the Courant-Snyder ellipses 1% (0) after 
the system passage are due to forth-order aberrations of the 
map A4 nc generated by the function S5. 

PROOF-OF-PRINCIPLE 
NUMERICAL EXAMPLE 

To complete the third-order apochromat design we need 
a drift-quadrupole cell with the coinciding periodic and 
apochromatic Twiss parameters and with the phase ad- 
vances satisfying the conditions ( fT~4-b and ( TTBT l. What is 
the minimum number of quadrupoles which are required 
for such a cell to exist? Using thin-lens approximation for 
quadrupole focusing one can show that this number can not 
be smaller than four, and with four quadrupoles we have 
many examples (found by direct numerical search using 
hard-edge quadrupole models) of the cells which satisfy all 
needed requirements. The parameters of one such cell are 
given in the Table 1 and its periodic (and apochromatic) 
betatron functions can be seen in Fig.l. The horizontal and 
vertical phase advances of this cell are equal to 144° and 
108° respectively, which means that the sequence of five 
such cells gives a third-order apochromatic beamline. 

Table 1 : Parameters of apochromatic four quadrupole cell. 
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The total number of quadrupoles in our third-order apoc- 
hromatic beamline is twenty, and it is interesting to com- 
pare its beam transfer properties with the beam transport 
through the chain of ten FODO cells which also has twenty 
quadrupoles and which is designed to have the same length 
and about the same average betatron functions. 

We chose for such comparison not one but three differ- 
ent FODO beamlines with the cell phase advances equal to 
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Figure 1 : Periodic betatron functions along apochromatic 
four quadrupole cell with yf x — 144° and [i c y = 108°. 

54°, 60° and 72° (see, for example Fig.2). The first and 
the third FODO beamlines have the phase advances which 
are equal to the half of the vertical and of the horizontal 
phase advances of our apochromatic four quadrupole cell 
respectively and, at the same time, are second-order apoc- 
hromats. The second FODO beamline is taken because, 
from one side, its phase advance 60° lies in between the 
phase advances of the two other beamlines, but, from the 
other side, its overall phase advance is not multiple of 180° 
and, therefore, it is not a second-order apochromat. 
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Figure 2: Betatron functions along two 54° FODO cells. 

Because it is not easy to find a single quantity which 
will characterize in a clear and general way the effect of 
aberrations which are nonlinear in the transverse variables, 
we will restrict our comparison to the study of the beam 
dynamics provided by the Hamiltonian 



H 



pi 



Pt 



+ (fc/2)-(x 2 -y 2 ) 



(16) 



after its linearization with respect to p 2 and p 2 i.e. we will 
look only at the effect of the energy deviation^ 

For each beamline, for each transverse plane and for 
each energy offset we track the nominal beamline Twiss 
parameters (Twiss parameters which are periodic for the 
nominal energy) through the system and, at the system exit, 



'if only chromatic aberrations are of concern and the influence of the 
geometric effects can be ignored, then the minimal number of cells re- 
quired for automatic cancellation of aberrations in our solution can be 
reduced from five to two because in this case only the resonances 
in 1 1 5t must be avoided. 
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Figure 3: Mismatch amplitudes for different beamlines. 

characterize the accumulated chromatic effects by calculat- 
ing first the mismatch parameter m p and then the mismatch 
amplitude 

a = m p + - 1, (17) 

which is a measure of the created optics distortion. The 
results of these calculations are presented in Fig. 3, where 
a x and a y denote the horizontal and vertical mismatch am- 
plitudes of the third-order apochromat and for each FODO 
chain only one curve is presented, because for the FODO 
beamlines the horizontal and vertical mismatch amplitudes 
coincide. One sees, as expected, that our third-order apoc- 
hromatic beamline demonstrates the best performance and 
the next comes the second-order apochromats with 54° and 
72° cell phase advances. 

Note that our theory is an asymptotic theory and uses 
Taylor expansion with respect to the phase space vari- 
ables. It certainly can fail at large energy deviations, that is 
demonstrated in Fig. 4. 
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Figure 4: Mismatch amplitudes for different beamlines. 
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